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The 2012 PDG reports a tension at the level of 3(7 between two exclusive determinations of 
\V u b\. They are obtained by combining the experimental branching ratios of B — s- TV and B nlv 
(respectively) with a theoretical computation of the hadronic matrix elements /b and the B — > n 
form factor f+{q 2 )- To understand the tension, improved precision and a careful analysis of the 
systematics involved are necessary. We report the results of the ALPHA collaboration for /b from 
the lattice with 2 flavors of 0(a) improved Wilson fermions. We employ HQET, including 1 jm\, 
corrections, with pion masses ranging down to i=s 190 MeV. Renormalization and matching were 
performed non-perturbatively, and three lattice spacings reaching cr 1 rj 4.1 GeV are used in the 
continuum extrapolation. We also present progress towards a computation of f+(q 2 ), to directly 
compare two independent exclusive determinations of \V u t,\ with each other and with inclusive 
determinations. Additionally, we report on preliminary results for /b s , needed for the analysis of 
B s -^pt+pL-. 
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1. Motivation 

The precise determination of the CKM matrix elements is a key for testing the Standard Model. 
Violations of CKM unitarity or discrepancies between independent determinations of the same 
matrix element can provide hints of New Physics. At the time when we started our work, a tension 
at the level of 3a between two exclusive determinations of | V u b| existed, as reported e.g. in the PDG 
of 2012. These determinations use the branching ratio (BR) for the processes B — > nlv and B — > TV 
from experiment combined with the form factor f+(q 2 ) and the B decay constant fy, respectively, 
from the lattice. Also an inclusive determination, based on a perturbative expansion in a s and an 
expansion in l/mj,, is possible [1]. The results reported by the PDG 2012, computed before ICHEP 
2012, can be summarized as follows [2]: 

| Kb | = 0.00323(31) (B->liv), |Kb| =0.00510(47) (B tv) , (1.1) 
| Kb | = 0.00441(34) (inclusive). 

At ICHEP 2012 the Belle collaboration reported a new result for BR(B — > tv) [3] based on a new 
set of data, obtained with a more sophisticated tagging of the B. This result, taken alone, would 
yield a value for |Kb| that is consistent with the exclusive determination from B — > n. However, 
more data and a careful inspection of the systematics involved are needed to draw more definitive 
conclusions. 

While the experimental precision in the differential decay rate for B — > %lv has by now reached 
good precision, B — > TV events are more difficult to reconstruct and there is an error of the order 
of 20% on the branching ratio. The situation on the theoretical side is the opposite: the lattice 
computation of a form factor is more challenging than that of a decay constant. 
At this conference we have presented the results for the determination of /b by the ALPHA col- 
laboration which use fully non-perturbative renormalization and matching, and CLS configurations 
with two degenerate dynamical quarks in the sea. A parallel effort to determine f+(q 2 ) in the same 
setup is ongoing: we have presented the progress reached so far, and the precision that we expect 
to achieve, once the full non-perturbative renormalization and matching at order 1 / nib have been 
completed. The comparison of these two exclusive predictions, in which the relevant hadronic 
parameters have been computed in the same setup, will provide a test as free as possible from sys- 
tematics. 

Recently LHCb has presented the first evidence for B s — > fi + jit~ , with a decay rate compatible with 
the Standard Model [4]. We have presented the determination by the ALPHA collaboration of the 
B s decay constant, /b s , that enters in the theoretical prediction of this decay. 

2. HQET 

Our computations were performed on CLS configurations, which have two degenerate 0(a) 
improved Wilson quarks. The ensembles used in this work have pion masses m n in the range 
190 Me V < m n < 450 MeV at three lattice spacings a, namely a G {0.078, 0.065, 0.045}fm. All 
of them have a spatial extent L such that m n L > 4, so that volume effects are expected to be very 
small. 

Even for our finest lattice spacing the b quark cannot be simulated directly, given that ami, > 1 ■ 
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However, for the low energy processes we are interested in, the mi, scale can be integrated out. Our 
approach is to use HQET, which is an expansion of the QCD Lagrangian in powers of I /mi,. At 
leading order the b quark is static, i.e. the Lagrangian involves no space derivatives. If we include 
terms up to order 1/m^ the Lagrangian becomes: 

-^HQET (x) = Wh M D Vh 0) - 0>kin ^kin (*) ~ ft>spin ff spin (x) , (2.1) 
^kin (X) = Wh (.X) D 2 Yh (x) , ^spin {x) = \j7 h (x) O ■ B \j/ h ( X ) . 

The corresponding expansion of the time component Ao of the heavy-light current (at p = 0) is 

s s 

, 4 tat = Wi YoYs ¥h , = Wi YsY Vh- (2-2) 

Since in HQET the 0(1 jm^) terms appear only as insertions in correlation functions, HQET is 
renormalizable order by order in I /mi, because the static theory is. Once the HQET parameters 

ft),- G {'«bare ) -Z^ ET ,c2\ft>kin,ft'spin} nave rjeen determined using non-perturbative matching [5], 
the continuum limit can therefore be taken safely. 

3. Matching 

The matching was performed in the Schrodinger Functional scheme in a small volume L\ pa 0.4 
fm, where am\, <C 1, and relativistic b-quarks can be simulated [6]. The HQET parameters can then 
be fixed by imposing the matching conditions for suitable observables <J> 

<J> HQET (Li,M,a)=<I> QCD (L 1 ,M,0), 0> QCD (Li,M,0) = ]im<t> QCD (L h M,a) , (3.1) 

a-s-0 

and hence the ft), inherit their dependence on the heavy-quark mass M from QCD. In particular here 
M is the RGI mass of the b-quark [7]. Using finite-size scaling recursively, we can then take the 
step L\ — > L2 = 2L\, and finally connect with large volumes L„ > max(2 fm,4/m w ). 

Having performed this matching procedure for each of the lattice spacings used in our large- 
volume simulations, we know the corresponding Nf = 2 HQET parameters ft);(Li,M,a) non- for a 
range of values of M in the neighbourhood of the b quark mass. 

4. Results 

The HQET energies and matrix elements are extracted at large Euclidean time t separations. 
The effects of excited states are exponentially suppressed like ~ e -( E 2- E i)t^ where E2 and E\ are the 
energies of the first excited state and the ground state, respectively. To achieve a better suppression, 
we solve the Generalized Eigenvalue Problem (GEVP) [8] 

C(t)v n (t,to)=Xn(t,to)C(to)v n (t,to), t <t<2t , (4.1) 

for anN xN correlator matrix C(t) with N = 3. Each entry of the matrix corresponds to a different 
Gaussian smearing level of the light quark field in the B -meson interpolating quark bilinear. The 
corrections to the energies and matrix elements so obtained behave like oc exp {— (E^+i —E\)i) and 



.HQET _ y HQET 
A 0,R — 



A A 
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HMxPT:/ B (zb,mps,a)/GeV / B (z b , m PS , o)/GeV /B.(zb,mps,a)/GeV 




Figure 1: Left: HM^PT extrapolation of /b; centre/right: linear extrapolation of /b and /b s . The blue, 
red and green points correspond to ensembles at a = 0.075 fm, 0.065 fm and 0.048 fm, respectively. The 
fit formulae evaluated at each given lattice spacing are shown in colour. The black curve is the chiral and 
continuum extrapolation. It is a fit to all the points shown in the figure plus other at the same pion masses 
but obtained using a different discretization for the heavy quark action. 

exp{ — (Eff+i —Ei)to} x exp{ — (£2 — E\){t — to)}, respectively. The residual systematic errors 
are kept under control by requiring a sta t > 3a sys . 

In phenomenological predictions, we need to know the HQET parameters at the physical mass 
of the b quark (Oj{L\,M^,a). To this end we impose m B (Li,Mb,m e ^ p \a = 0) = m£ p = 5279. 5MeV. 
The mass of the B at physical pion mass «% xp is obtained through a chiral and continuum extrapo- 
lation [9]: 

m B (z,m K ,a,n)=B(z)+Cml-—^-,mi + D n a 2 , g = 0.51(2) [10]. (4.2) 

In HQET the B-meson mass is given by ra B = m bare + £ stat + a\mE kin + co sp i n £' spin , so that 

M b = 6.56(15) (06) z GeV, or equivalently mjp(m b ) =4.22(10)(4) z GeV . (4.3) 

For the following analyses, we use the values of HQET parameters obtained from an interpolation 

to z = Zb- 



4.1 B and B s decay constants 

Using the HQET formulas to combine the matrix elements and energies with the matching 
parameters at z = Zb we can compute /b and /b s through a chiral and continuum extrapolation 
dictated by HMChPT [11, 12]: 

3 1+?" 2 



f B (m n ,a,n) =b 
Our analysis for ICHEP 2012 gives 



1 



4 (4nf„) 



^m n \n{m n ) 



9 9 

+ cm n + d n a . 



(4.4) 



/ B = 193(9) stat (4) z MeV, f Bs = 219(12) stat MeV , 

where the error coming from the chiral extrapolation is determined by comparing to a linear ex- 
trapolation in m n (see Fig. 1). For /b s not all ensembles are analysed yet. For more details, see 
[13]. Our value for /b is compatible with the values found by other collaborations [14, 15, 16, 17]. 



4 



|y„f,| in lattice QCD 



F. Bernardoni 



4.2 B — > K form factor 



The form factor f+(q 2 ) is defined through the Lorentz decomposition of the matrix element: 



(K{p n )\V»\B(p B ))=f + {q 2 



2 _ 2 



Pb+Pk 



+ fo(q 2 ) 



mi — m 



(4.5) 



where = p 1 ^ — p^- To extract this matrix element from lattice simulations, we consider the ratio: 



I e- ipx * (Pi, (t„ + tB,x 7[ )V» (f B , xb )P M (0) ) 
I e-fr** (Pa (x„)P u (0) ) I (P h i (x B )P U (0)) 



<^(p ff )|V|fl(pB)) 



hm /?(^,fB)^ /2+mB ' B/2 



(4.6) 



where Z 5 // and P^ are interpolating operators for the % and the 6 meson, respectively. 
One additional difficulty in comparison to 
the extraction of decay constants is the pres- 
ence of large finite T effects, where T is 
the temporal extension of the lattice. So far 
we have restricted to few lattices and to the 
static limit to demonstrate the feasibility of 
this computation in our setup. Our results 
show that finite T effects can be understood 
in the transfer matrix formalism (see Fig. 2) 
and that a precision of 5% is achievable. 
Combining our result at largest q 2 , q = 
Pti — Pb, with the experimental data available 
and a parametrization of the q 2 dependence 
based on very general properties like uni- 
tarity and analyticity [18], \V u b\ is obtained 
with a 15% precision. This does not include 
the necessary extrapolation in the light quark 
mass and the lattice spacing. These missing 
steps will soon be coming out. For more de- 
tails, see [19]. 




Figure 2: t K dependence for the considered ratio 
eq. (4.6) on one of our lattices with m K w 300 MeV. 
We use a smeared interpolating operator for the B and 
the K. The red curve is obtained by fitting the finite 
T effects with the prediction from transfer matrix the- 



ory: R(t K ,t B ) 
rani etc is 



-E n (T-1t n ) 



y/l +e -E n (T-2tx) 



with A and B as pa- 
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